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求解 多 目标 带 时 间 窗 VRP 问题 的 文化 狼 群 算法 


李 人 小川 ， 刘 媛 华 ， 王 影 歌 
(上 海 理工 大 学 管理 学 院 ， 上 海 200093) 


摘 要 : 针对 以 最 小 化 总 距离 和 车 辆 数 为 目标 的 带 时 间 窗 车 辆 路 径 问题 的 求解 , 提出 一 种 基于 文化 基因 的 狼 群 算法 。 
根据 客户 近邻 度 和 矩阵 构建 初始 狼 群 ， 提 高 狼 群 的 优良 性 。 结 合 问题 特征 与 狼 群 算法 的 寻 优 思想 ， 重 新 定义 了 群体 的 
游 走 和 突袭 行为 。 提 出 头 儿 变异 策略 ， 增 加 了 最 优 解 的 令 域 搜索 范围 。 受 文化 算法 启发 ， 将 当前 最 优 解 的 有 利信 息 
作为 信仰 空间 的 知识 来 指导 狼 群 进化 方向 ， 提 高 搜索 效率 。 通 过 仿真 实验 及 与 其 他 算法 对 比 ， 表 明文 化 狼 群 算法 求 
解 车 辆 数 更 少 、 总 距离 更 短 ， 且 具有 更 好 的 稳定 性 。 
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Cultural wolf pack algorithm for solving multi-objective VRP with time window 


Li Xiaochuan, Liu Yuanhua, Wang Yingge 
(Business School, University of Shanghai for Science & Technology, Shanghai 200093, China) 


Abstract: To solve the vehicle routing problem with time window for the shortest length of vehicle travel and the minimum 
number of the using vehicles, this paper proposed a cultural wolf pack algorithm. To improve the quality of the wolf pack, 
this paper constructed the initial wolf pack according to the customer's nearest neighbor matrix. Considering the features of 
this problem and the wolf pack algorithm, this paper redefined the intelligent behaviors of artificial wolves. The mutation of 
leader wolf increased the searching range of the optimal solution. According to the favorable information of the current 
optimal solution, this paper constructed the knowledge of the belief space to guide the evolution direction of the wolf pack 
and improve the search efficiency. Simulation results show that the cultural wolf pack algorithm can get less number of 
vehicles, shorter total distance and better stability than other algorithms. 

Key words: vehicle routing problem; time window; multi-objective; cultural wolf pack algorithm 


0 引言 吸取 经 验 来 指导 自身 或 他 人 行为 的 启发 而 提出 ， 对 智能 算法 
喇 的 自 组 织 过 程 赋予 显 性 指导 条 件 。CA 为 上 层 信仰 空间 (也 
车 辆 路 径 问 题 (vehicle routing problem ,VRP) 具有 非 线 称 信 度 空间 ) 和 底层 种 群 空 间 组 成 的 双 层 并 行 结构 59]， 信 仰 
性 、 非 凸 性 、 复 杂 性 和 约束 条 件 多 且 相 互 之 间 难 以 协调 的 特 空间 以 接受 函数 的 形式 记录 种 群 迭代 过 程 中 的 有 用 信息 ， 并 
点 ， 是 典型 的 非 线 性 规划 模型 ， 也 被 证 明 为 NP-hard 问题 。 对 其 分 析 总 结 形成 知识 来 指导 种 群 空 间 中 粒子 的 后 续 进 化 过 
该 问题 自 1959 年 由 Dantzig 等 人 后 提出 以 来 一 直 受 到 众多 学 程 。CA 为 智能 算法 的 改进 提供 了 一 种 思想 ， 各 种 符合 CA 
CS 者 的 研究 和 关注 ， 拥 有 更 多 限制 条 件 的 车 辆 路 径 问 题 也 相继 结构 的 算法 都 可 以 作为 种 群 空间 的 一 部 分 参与 进化 0， 从 而 
被 提出 。 带 时 间 窗 车 辆 路 径 问 题 (VRP with time window, 是 高 基本 算法 的 性 能 。 狼 群 算法 (wolf pack algorithm, WPA ) 
VRPTW) 由 SavelsberghD 提 出 ， 其 在 VRP 的 基础 上 增加 了 模拟 大 自然 中 狼 群 分 工 协作 、 和 群体 捕食 的 行为 进行 迭代 寻 优 ， 
客户 接受 服务 时 间 区 间 的 约束 ， 主 要 优化 目标 为 最 小 化 车 辆 。” 符合 文化 算法 的 基本 结构 中 。 本 文 提出 一 种 文化 狼 群 算法 
数目 、 最 小 化 总 行驶 距离 、 最 小 化 总 成 本 等 B]。 (cultural wolf pack algorithm，CWPA ) 用 于 求解 多 目标 
VRPTW 比 VRP 更 符合 现实 的 物流 配送 过 程 , 但 其 求解  VRPTW, 用 狼 群 算法 来 进行 迭代 寻 优 , 提取 每 代 最 优 路 径 ( 头 
过 程 也 更 加 复杂 ， 智 能 启发 式 算法 在 求解 该 问题 时 具有 一 定 ， 狼 ) 中 的 最 优 车 辆 基因 段 作为 形势 知识 ， 并 以 其 子 段 作为 规 
的 优越 性 ， 受 到 很 多 学 者 的 青睐 。 目 前 已 有 多 种 智能 优化 算 ” 范 知 识 来 指导 探 狼 游 走 和 猛 狼 突袭 行为 , 避免 进化 的 盲目 性 
法 应 用 于 VRPTW 的 求解 , 如 遗传 算法 差分 进化 算法 个 、 ”从 而 提高 搜索 效率 。 通过 实验 验证 了 本 文 CWPA 较 蚁 群 算法 
蚁 群 算法 [由 等 。 以 上 算法 均 对 VRPTW 问题 的 求解 进行 了 蝙蝠 算法 、 禁 忌 搜索 算法 等 对 比 算法 在 求解 多 目标 VRPTW 
一 定 程度 的 探究 ， 但 是 遗传 算法 和 差分 进化 算法 存在 易 早 熟 问题 时 的 优越 性 ， 为 该 问题 的 求解 提供 了 一 种 新 思路 。 
收敛 的 问题 ， 蚁 群 算法 易 陷 入 局 部 循环 ， 人 迄今 为 止 尚 没 有 哪 _ i 
种 算法 对 VRPTW 问题 具有 绝对 的 求解 优势 ， 因 此 寻找 对 该 1 ”多 目标 VRPTW 数学 规划 模型 


一 


TE 


a 


问题 求解 性 能 更 优 的 算法 仍 具 有 重要 的 现实 意义 。 1.1 问题 描述 
文化 算法 (cultural algorithm, CA ) 由 Reynolds023 受 人 类 本 文 研 究 的 多 目标 VRPTW 问题 可 以 描述 为 由 k 辆 规格 
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发 ， 服 务 完 该 路 径 上 


李 小 川 ， 等 


等 : 求解 多 目标 带 时 间 窗 VRP 问题 的 文化 狼 群 算法 


围攻 捕食 等 行为 ， 根 据 狼 群 的 食物 分 配 规则 来 进 
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行 种 群 更 新 。 算 法 的 初始 释放 狼 群 进行 搜索 ， 选 取 嗅 到 食物 


路 径 ， 使 得 配送 车 辆 


相同 的 车 为 (n>#) 个 客户 配送 货物 ， 每 辆 车 从 配送 中 心 出 。” 狐 狼 突袭、 
的 客户 后 回 到 配送 中 心 ， 通 过 合理 规划 
数 和 总 距离 最 小 。 同 时 还 要 满足 以 下 条 

牛 : a) 每 辆 车 装载 的 总 货物 都 必须 满足 其 最 大 载重 量 限制 ; 


b) 每 辆 车 可 以 服务 多 个 客户 , 但 每 个 客 


i 只 在 时 


分 别 为 客户 i 接受 服务 的 最 


户 都 只 能 
| 间 区 间 [Zi,7] 内 接受 服务 ， 划 
早 和 最 晚 时 


司 ， 若 配送 车 辆 早 于 


T; 到 达 ， 则 需要 等 到 7 进行 服务 ,但 最 大 等 待 时 间 要 受到 约 


或 者 返回 配送 中 心 。 
1.2 参数 定义 


N= 亿 2… 吕 :所 有 客户 的 集合 


K ={1,2,…, 


} :所 有 配送 车 辆 的 外 


ly 


区 


必 : 客 户 ; 与 客户 /之 间距 离 


W :车 辆 的 最 大 载重 量 


i :客户 i 的 货物 


:客户 i 的 


需求 量 


页 大 服务 的 客户 点 集 


时 间 窗 


7 :车 辆 到 达 客户 i 处 的 时 间 点 
7 : 辆 在 客户 ; 处 的 等 待 时 间 


Si :客户 的 服务 时 间 


决策 变量 : 


1 车 辆 K 
Xi = 
0，others 


1, 
Xix = 0 


1.3 数学 模型 


R 务 完 客 户 ;后 接着 服务 7/ ; 


束 ; d) 车 辆 k 在 客户 i 处 的 服务 结束 后 即刻 去 往 下 一 客户 处 


jeN (1) 
客户 由 车 辆 k 服 务 
ieN,keKk 
， others » "© 。 O) 
NK 
CI =min》, 》 区 (3) 
1 ka 
上 : 丰 蜂 
C, =inin DY 3 xdy (4) 
i=l j= k=l 
NK 
2 -1, jeN (5) 
i=] k=l 
N 
Drixx SW ,kek (6) 
i=l 
天 
> xx =1,ieN (7) 
k=1 
NV N 
的 = pe ,keKk (8) 
问 半 
N N 
De ,keKk (9) 
i=2 j=2 
x =|P|-1 
站 不 I2| ,KE 天 (10) 
Tu<TitTt+S STi,ieN (11) 
W =max{0,T,—T,} ,ieN (12) 


式 (3) 和 “4) 分 别 表示 车 辆 数 最 少 和 总 距离 最 小 的 目 


标 函 数 ; 式 (5) 约 束 了 每 辆 车 
表示 车 辆 装载 货物 不 能 超过 该 车 最 大 载重 


气味 浓度 最 大 〈 即 适应 度 值 最 好 ， 


取 值 fs ) 的 作为 头 狠 ， 


另 选 除 头 狼 外 适应 度 值 较 好 的 W 只 狼 作为 探 狼 , 每 只 探 狼 有 


灵 次 游 


置 ， 


机 会 ， 即 在 个 不 同方 向 进行 局 部 搜索 ， 游 走 过 程 
中 若 嗅 到 的 食物 气味 浓度 /大 于 


Bm ， 则 该 探 狼 取 代 头 狼 位 


Fe = 三， 否则 退 


口 


lh 


原 位 ， 继 续 游 走 。 头 狼 发 出 召唤 


? 


猛 狠 向 头 狠 方 向 奔 克 ， 


奔袭 途中 嗅 到 的 食物 气味 浓度 若 大 于 5， 则 


为 头 狼 并 发 起 召唤 ， 否 


接着 狐 旬 


度 值 
淘汰 。 


大 到 小 排列 ， 


狼 群 的 游 走 、 


A 


走 过 程 中 存在 


于 / 司 


搜索 能 力 ， 


良 联 同 探 狼 对 猎物 发 起 车 
到 的 食物 气味 浓度 的 大 小 决定 是 否 更 新 位 置 。 将 狼 群 按 适 应 
选取 前 $ 只 


奔袭 行 
4 用 攻 行 为 进 
WPA 在 寻 优 过 程 中 主要 存在 以 下 两 方 


此 过 程 类 似 粒 子 群 算法 的 社会 行为 ， 
该 猛 狼 立即 转 
则 继续 奔袭 至 与 头 狼 距离 4<mind 。 
攻 ， 根 据 其 围攻 结束 后 若 嗅 


果 留 到 下 一 代 ， 其 余 的 则 被 


为 使 得 狼 群 算法 拥有 较 好 的 局 部 和 
步 增 强 了 算法 寻 优 性 能 ， 但 是 
甸 不 足 :，a) 探 狼 在 游 


定 的 盲目 性 ， 


学 习 能 力 较 差 ; b) 狐 狼 在 奔袭 


时 的 距离 判别 法 存在 较 大 的 主观 性 ， 且 增加 了 算法 复杂 度 
[1617,6CA 的 重要 思想 为 将 种 群 进化 过 程 中 的 经 验 提 取出 来 ， 


经 过 分 析 处 理 形成 知识 ， 并 存储 于 信仰 空间 。 


信仰 空间 中 的 


知识 一 般 超 越 算 法 的 自学 习 能 力 ， 从 而 使 算法 具有 跟 更 高 的 
E。 本 章 提出 的 CWPA 首先 由 WPA 进行 路 径 搜 索 ， 将 最 


智 色 


优 路 径 中 的 最 
狐 狼 的 行为 进 


2.1 


眠 


包车 辆 信息 作为 信仰 空间 中 的 知识 ， 对 探 狼 和 
行 指导 ， 增 加 狼 群 的 寻 优 性 能 。 
适应 度 函 数 的 构造 


采取 精英 保留 策略 ， 以 路 径 总 距离 和 车 辆 数 作为 评价 指 
标 ， 构 造 如 下 适应 度 函 数 : 


™ 


数 ; 


2.2 


Pp: DD; 为 第 i 个 解 的 路 径 总 距离 ; 
4 的 加 入 是 为 避免 极端 解 〈 比 如 车 辆 数 
较 多 ， 但 总 距离 较 短 的 解 ) 的 影响 。 


4 为 权重 系数 ， 


算法 编码 


f 


G 
C, 


mE D 于 
5D (13) 


ieN 


GC; 为 第 i 个 解 所 需 车 辆 


CWPA 采用 整数 编码 来 表示 车 辆 路 径 ， 配 送 中 心 


j 0 表 


不 ， 


客户 则 ) 


AEt 


j 1 从 配送 


车 有 


2.3 


次 只 能 服务 一 个 客户 ; 式 (6) 


三 | 
里 ; 


式 (7) 约束 了 


每 个 客户 只 能 被 一 辆 车 服务 ; 式 (8) 表示 每 辆 车 从 配送 中 心 


出 发 ， 最 后 还 要 回 到 
程 的 连续 性 ， 即 服务 


同 约束 了 客户 只 在 时 间 窗 内 接受 服务 。 


2 求解 多 目标 VRP 
狼 群 算法 将 狼 群 的 扣 


配送 中 心 ; 式 (9) 约束 了 车 辆 上 配送 过 
完 客户 i 后 也 必须 从 该 点 出 发 ; 式 (10) 
约束 了 每 辆 车 的 路 径 都 不 存在 子 回 路 ; 式 (9) 和 式 (10) 共 


TW 的 文化 狼 群 算法 
食 过 程 抽象 为 探 狼 游 走 、 头 狼 召 唤 、 


初始 狼 群 的 构造 
1) 路 径 选 择 规则 


个 客户 进行 配送 , 序列 03501420 表示 
FP 心 出 发 ,依次 服务 客户 3 和 5 后 回 到 
有 2 从 配送 中 心 出 发 , 依次 服务 客户 1、4 和 2 后 下 
中 心 ， 基 因 位 4 中 的 0 表示 车 辆 1 的 终点 和 车 辆 2 


其 对 应 的 编号 表示 。 例 如 由 1 个 配送 中 心 为 5 


两 辆 车 进行 配送 , 车 
配送 中 心 ， 
到 配送 
的 起 点 。 


< 


初始 狼 群 个 体 对 路 径 节点 的 合理 选择 ， 可 以 提高 狼 群 的 
整体 质量 ， 加 快 收敛 速度 。 设 di 为 客户 i 与 i 之 间 的 距离 ， 


且 =d;， 由 此 构 如 下 近 领 度 和 矩阵 : 

COD)=(00)w ,isjeN (14) 
其 中 : 0y e {2,…， nn 由 ,表示 客户 i 与 客户 j 的 近邻 度 ( ci=% )， 
如 cxs =5 表示 客户 2 是 距 客 户 6 第 5 近 的 客户 。 


若 车 辆 当前 所 在 客户 节点 为 i， 选 取 客户 j 作为 下 一 节 


点 概率 PG,j) 可 由 式 (15) 得 到 
el EA 
270,)) 
PD)=1 症 , ,bjeN (15) 
(人 > Cnm 
六 nD) 
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(16) 


jeN, jzi Ci 
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内 中 7(,j)= 
6 为 最 大 近邻 度 关 

度 小 于 一 定 范围 

表示 与 客 


7 j) 


“人 mij 被 选 为 下 一 有 


i 的 近 
务 点 。 当 
个 客户 都 上 


2) 初始 解构 造 步骤 
设 配 送 车 辆 规格 相 


载重 量 ， 


a) 初始 化 k=1， 


客户 i 被 选 为 


R 务 点 。 为 了 保持 种 群 多 样 


定常 数 , 在 车 辆 路 径 选 择 中 认为 近邻 
内 的 客户 作为 相 邻 节点 较为 合理 ; 
户 的 近邻 度 6 的 客户 了 将 以 较 大 的 概率 


四 =0.98 ， 


性 ， 与 客户 


邻 度 >co 的 客户 也 


只 被 服务 一 次 。 


服务 点 后 ， 即 刻 令 Pq.D=0， 


同 ,载重 量 为 W ，w 为 车 辆 的 剩余 


了 机 会 以 较 小 概率 被 选 为 下 一 服 
以 保证 每 


到 为 车 辆 上 服务 的 客户 点 集合 ，K 为 最 多 使 用 的 车 


辆 数 ， 车 辆 从 配送 中 心 


w=W, I={0}。 
b) 随机 产生 车 辆 1 的 初始 服务 客户 5 中 w=w 

L=LU{s}, 邻 Pi,s)=0。 
c) 生成 一 个 [0 之 间 

辆 的 配送 节点 i。 


的 随机 数 P， 


出 发 时 间 无 先后 顺序 ， 则 初始 解 的 构 


74， 
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原则 上 ， 配 送 车 量 在 不 超载 的 
好 ， 如 果 装 载 货物 过 少 则 存在 较 大 
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前 提 下 装载 的 货物 越 多 越 


的 浪费 ， 


与 行驶 距离 两 方面 的 因素 确定 最 优 车 辆 的 判定 准则 如 式 (17) 


所 示 ， 使 sow,Do 取 最 小 值 时 的 车 辆 即 为 最 优 车 辆 。 


g (Wi, D = 
其 中 : wi 为 车 辆 k 的 剩余 装载 量 ; 
2) 信仰 空间 的 更 新 
信仰 空间 的 更 新 采取 精英 保留 
[XQ©O— X,Y K,(t+l) 


| 


其 中 : X(t+lD) 为 第 1+1 代 的 


势 知 


优 车 辆 服务 客户 集合 ; 
gi 为 前 1 代 的 最 小 8 值 ，K.Q 为 
2.5 ” 探 狼 游 走 


xD. 


(17) 


Di 为 车 辆 的 行驶 总 距离 。 


策略 ， 


XQ 中 所 有 客 


选取 除 头 狼 外 适应 度 值 较 好 的 M 只 人 工 狼 


Qa+l o 


狼 比 例 因 子 。 探 狼 的 游 走 
其中 PQ 为 第 : 代 规 范 知识 中 的 基 


m 取 | 立 | 内 的 随机 


根据 式 〈15) 选择 和 


让 


d) 若 q<Sw, HTitTi sh, Nw=w-q, L=Luf, 


令 PG,i)=0， 
e) k= 
日 将 
g) 重复 执行 
h) 中 4| 为 集合 上 

示 车 辆 k 完成 配送 任务 


转 c)， 
k+l, 
L. 整理 


b)~ 们 ? 


否则 执行 e)。 
若 k>K， 则 令 k=K， 转 c)。 
为 可 行 解 


加 | 王 Ilo] 。 
成 YX 个 初始 解 。 


素 的 顺序 排列 ， 序 列 末尾 的 0 表 


a 
了 妆 | | 必 
再 | Ft [区 


到 配送 中 心 。 


2.4 信仰 空间 
1) 信仰 空间 的 形成 
s 间 对 每 次 迭代 中 最 优 个 体 的 优秀 基因 进行 分 析 构 


信仰 空 


成 形势 空 


[19] 。 


车 辆 路 径 


zs 间 和 规范 空间 ， 为 种 群 的 进化 方向 及 范围 提供 指导 


问题 的 可 行 解 由 每 辆 车 服务 的 客户 编号 排列 组 


成 ， 因 此 可 以 提取 初始 种 群 最 优 个 体 中 表示 最 优 车 辆 配送 的 


客户 集合 
知识 。 


例如 


作为 


程 分 别 如 


|5,3| 三 个 基因 段 。 


， 则 形 
图 1(a) 和 (b) 所 示 , 其 中 规范 知识 包括 |2,5,3|、 


攻势 知识 ， 其 部 分 连续 配送 客户 子 集 作 为 规范 
初始 种 群 最 优 个 体 的 最 优 车 辆 (以 下 简称 最 优 车 
辆 ) 配送 客户 集 为 {[2.5, 引 


势 知 识 和 规范 知识 的 形成 过 
2,5| 和 


形势 知识 
3 5 5 
(a) 形 势 知识 形成 过 程 


(a) Formative process of Situational knowledge 


规范 知识 


y vl 
Pp 3 3 7 
全 


(b) 规 范 知识 


成 过 程 


(b) Formative process of normative knowledge 


图 1 信仰 空间 


成 过 程 


Fig.1 Formatiive process of belief space 


次 数 根 提 


整数 ，W 为 狼 群 个 体 总 数 ， 


昌 和 迭代 次 数 动态 取 值 ，，， 
段 个 数 。 探 狼 游 走 


因 


中 ， 首 先 将 探 狼 基因 序列 中 的 “0” 移 除 ; 然后 将 规范 知识 中 


攻 势 知识 更 新 如 下 : 
，So(t+D< gD) (18) 
XU K+) , gs(t+D) > 8 & KO MK, t+h)=G 

识 ; X,G 为 第 
识 中 8 值 最 大 的 车 辆 服务 的 客户 集合 ;Ks(t+tD) 为 第 t+1 代 最 
sodt+D 为 第 1+1 代 最 优 


1 代 形 势 知 


车 辆 的 8 值 ; 


户 点 集合 。 


作为 探 狠 ， 
2 为 控 


过 程 


的 每 个 基 医 


段 分 别 在 其 首尾 和 中 间 


随机 位 


各 添加 


次 ， 每 


次 添加 需 将 原 序列 中 的 相应 编号 移 除 ， 例 如 某 探 狼 的 除 “0” 


基 医 


序列 为 52.8.9.47 ， 


在 其 末尾 添加 8,4， 


最 后 将 其 转换 为 可 行 解 编码 。 若 该 探 狼 嗅 到 食物 


eo 则 取代 头 狼 位 置 ， 
续 游 走 至 Ab 次 。 
2.6 ” 猛 狼 突 效 

特 淮 辣 尖 积 力 同 突 大 让 二 下 

围 ， 提 高 了 搜索 效率 。 以 其 中 一 只 


并 发 起 召唤 ， 否 则 退 


加 了 当前 最 
狐 狼 为 例 ， 


多 解 的 领 
具体 步骤 为 : 


得 到 8,4,5,2,9,7 ; 


气味 浓度 大 
回 原 位 ， 继 


i 域 范 


a) 令 k=1, 提取 规范 知识 中 表示 8 值 最 小 车 辆 所 有 配送 


客户 的 序列 L。 
b) 根据 近邻 度 和 矩阵 分 别 找到 
近 的 客户 ab 。 
o) 将 蕊 插入 到 猛 狼 序列 Mw 


与 必 首 


d) 移 除 MW' 中 与 5 重复 的 客户 得 到 


尾 客户 a,b 距离 最 


ha 之前， 得 到 序列 MW'。 


序列 MW? 。 


e) 将 Mr 中 的 请 移 到 2 之 后 得 到 猛 狼 完成 突袭 后 的 序 


列 MW3 。 
f 计算 适应 度 f。 。 
g) 若 刻 ,小 于 当前 头 狼 适 应 度 方 ， 执 行 bj; 和 否则， 退回 
原 位 后 转让 )。 
h) 该 独 狼 取代 头 狼 ， 更 新 信仰 空间 ， 并 重新 发 起 召唤 。 
iD 令 k=k+1， 若 k>2 ， 该 狼 结 束 奔袭 ， 若 上 >2 ， 令 
c=b',bi=al,al=c ， 转 Cc)。 
2.7 头 狼 变异 
CWPA 用 头 狼 变异 行为 取代 WPA 中 的 狼 群 围攻 行为 。 


WPA 中 群 狼 围攻 行为 实质 上 是 各 粒子 对 当前 最 优 粒子 的 邻 


进行 狼 群 围攻 


己 经 充分 利用 了 头 狼 信息 ， 如 果 再 
许多 无 效 操作 。 另 外 ， 如 果 将 每 只 
次 变异 操作 ， 则 狼 群 


际 上 
所 以 会 节省 算法 运行 时 间 。 头 狼 变 
异 对 象 为 单个 客户 点 ， 例 


异 采 取 单 点 变异 策 
如 头 狼 的 去 “0” 序 列 为 


域 进行 深度 搜索 的 过 程 , 而 本 文 CWPA 的 改进 狐 狼 突袭 操作 


操作 会 增加 


工 狼 的 围攻 行为 看 做 一 
围攻 相当 于 进行 了 N1 次 变异 
3.1 节 可 知 头 狼 变 异 次 数 为 3 时 即 可 取得 较 好 的 结果 ， 


行为 ， 实 


已 和 


各 ， 变 


2,5,9,3,6,8,4,1,7 ， 


客户 2,5 由 车 辆 1 配送 ， 客户 93,6,8 | 


车 辆 2 


配送 ， 客 户 417 由 车 辆 3 配送 ， 图 2 表示 头 狼 变异 过 程 。 
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2.8 CWPA 算法 步骤 


狠 群 适应 度 ， 取 最 


a) 初始 化 各 参数 ， 设 置 初始 人 了 
b) 根据 3.3 节 中 步骤 进 


图 


2 头 狼 变 异 过 程 


Fig.2 Mutation process of leader wolf 


[ 狼 个 数 。 
行 狼 群 初始 化 ， 
优 适 应 度 值 的 人 工 狼 被 选 为 头 狼 。 


由 式 (13) 计算 


c) 根据 3.4 节 内 容 构 造 或 更 
[NN 
[latl'o 


d) 随机 选 一 


新 信仰 空间 。 
| 内 的 整数 作为 探 狠 个 数 ， 计 算 第 


代 每 只 探 狼 的 游 走 次 数 Ab 。 若 游 走 过 程 f( 探 狼 i) < hh ， 该 探 


转 


2.9 ”CWPA 算法 特点 分 析 
是 出 的 CWPA 了 
探 狼 个 体 在 各 自 的 邻 域 范围 
无 解 附近 , 而 全 局 最 优 
良 个 体 


域 最 
狐 狼 突袭 阶段 ， 狸 间 
前 头 旬 
的 搜索 邻 


的 


狼 取代 头 狠 位 
走 至 xb 次 。 
据 3.6 节 步 又 进行 独 狼 突袭 。 
人 执行 头 狼 变异 操作 ， 
良 取代 头 狼 位 置 ， 转 g)。 
g) 更 新 狼 群 , 记录 当前 头 狼 适 应 度 值 。 若 1<G , 则 :=t+l， 


则 该 变异 


e) 根 


c); 否则 算法 结束 ， 输 出 最 优 结果 


本 文 


EL ， 转 c) 更 新 信仰 空间 ， 再 转 e)， 否 则 继续 游 


若菜 次 变异 后 适应 度 值 优 于 有 ， 


o 


要 分 为 三 个 阶段 : 探 狼 游 走 阶段 ， 


或 ， 也 有 
搜索 ， 狼 群 基本 


内 进行 深度 
解 也 有 较 大 概率 在 这 些 点 中 产生 ; 


局 部 搜索 ， 处 在 J 


局 


不 


集中 在 最 


利用 了 文化 空 


间 中 的 有 利信 息 向 当 


良 靠近 ， 避 免 了 突袭 的 盲目 性 ， 相 当 于 


曾 加 了 最 优 个 体 


利于 种 群 的 全 局 进化 方向 ， 经 过 前 两 阶段 
优点 附近 ， 这 时 执行 头 狼 变异 行 


有 利 了 


首 加 种 群 多 样 性 。 所 以 CWPA 在 


为 
局 部 和 
3 


为 了 验 记 
VRPTW 的 效能 ， 利 月 


加 


局 寻 优 性 能 。 
仿真 实验 结果 及 分 析 
FE 本 文 所 提 文 化 狼 群 算法 CWPA 求解 多 目 


里 论 上 具有 较 好 的 


标 


MATLAB R2012a 软件 进行 两 组 仿真 


实验 ， 仿 真 环境 为 处 理 器 Intel® CoreGM i5-2350M CPU 


2.30 GHz,RAM 4 GB,Win7 
CWPA 算法 参数 设置 
CWPA 与 相 
性 原则 , 旨 在 使 
针对 CWPA 部 分 导 
集 上 进行 了 多 组 


3.1 


c 为 5~7， 探 狼 比 作 


范 知识 


区 


续 客 


户 子 人 自 


形成 机 制 
3 所 示 的 规律 。 医 
配送 的 客户 数 ，S-1 
辆 配送 客户 子 集中 


关 文 献 中 


™ 


操作 系统 的 PC 机 。 


他 改进 算法 的 公共 参数 采取 一 致 


大 | 


其 运行 结果 与 文献 中 的 实验 结果 更 具 可 上 


要 的 特有 参数 , 本 文 在 Solomon 标准 
测试 实验 ， 可 以 得 出 : 最 大 近邻 度 判定 常数 
| 因子 x 为 4 或 5， 头 狼 变异 次 数 
时 能 取得 相对 较 好 的 求解 效果 。 此 外 ，CWPA 对 信仰 
的 灵敏 度 也 较 大 ， 分 析 多 组 实验 结果 


1 为 3 
空间 
总 结 出 


区 


中 8 表示 经 式 (17) 判定 的 最 
，S-2，…2 表示 规范 知识 中 提取 的 最 优 车 
包含 的 客户 数 ， 取 8 时 表示 不 执行 规范 知 
识 操 作 ， 取 5S-1 时 表示 规范 知识 提取 包含 S-1 个 客 
发 ， 取 5S-2 时 表示 规范 知识 


车 辆 所 


户 数 的 连 
提取 包含 S-1 和 5-2 个 


客户 数 的 连续 客户 子 集 ， 同 理 取 2 时 表示 规范 知识 提取 包含 


S-1、S- 


搜索 效 


2、 


明 
个 


…2 个 客户 数 的 连续 客户 子 集 。 可 以 看 出 
较 不 执行 规范 知识 时 


取 S-1 时 
取 S-2 或 更 小 


a 


明显 提 


值 时 对 搜索 效果 
呈 非 线性 增加 ， 


大 | 


的 连续 客户 子 集 。 


的 贡献 度 基本 为 零 ， 但 是 会 使 算法 运行 时 间 
此 本 文 规范 知识 只 


提取 包含 S-1 个 客户 数 


图 


3 规范 知识 运行 机 制 对 算法 性 能 影响 


Fig.3 Relation between operating mechanism of normative 


knowledge and algorithmic performance 


3.2 仿真 实验 1 


下 


L 送 中 心 为 20 个 客 


首先 以 文献 [3] 中 的 算 例 为 例 进行 
法 、 文化 狼 群 算法 和 文献 [3] 中 提出 的 单 点 重组 精英 遗传 混合 
蝙蝠 算法 MRGBA 与 文献 [11] 中 所 提 的 单亲 遗传 混合 蚁 群 
法 PHACA 求 得 的 结果 进行 对 比 。 算 例 可 
户 配 送 货物 ， 配 送 车 玫 
重量 均 为 8 ,车速 均 为 40km/h, 车辆 在 客户 点 的 最 大 等 待 时 


仿真 实验 ， 并 对 狼 群 算 


= 


以 描述 为 : 1 


外 
个 


规格 相同 , 最 大 载 


间 为 45h ， 求 解 目标 为 最 小 化 配送 车 辆 和 


如 表 1 所 示 〈 其 中 客户 0 为 配送 中 心 )。 


总 距离 ， 客 


户 信息 


表 1 客户 信息 表 
Table 1 Clients information 
客户 ”客户 坐标 ”需求 量 最 早 接受 最 晚 接受 服务 
编号 /km At 服务 时 刻 几 服务 时 刻 必 ”时间 笑 
0 (70,70) 0 0 0 0 
1 (107,77) 3.4 1.5 3.5 0.2 
2 (109,139) © 0.8 3.0 $5.5 0.2 
3 (120,22) 3.9 2.0 5.0 0.2 
4 (48,47) 1.9 1.5 5.0 0.4 
5 (116,22) 3.2 2.0 5.5 0.2 
6 (12,138) 1.4 4.0 8.0 0.5 
这 (86,40) 22 1:3 3.0 0.4 
8 (121,124) 2.1 2.5 5.0 0.2 
9 (61,57) 3.5 4.5 8.0 0.2 
10 (40,113) 2.3 2.0 3.5 0.2 
11 (129,24) 1.8 2.5 $5 0.4 
12 (12,84) © 1.6 2.0 4.5 0.4 
13 (44,116) 2;7 3.0 6.0 0.2 
14 (102,52) 1.5 3.5 5.5 0.2 
15 (41,36) -1.3 4.0 7.0 0.2 
16 (132,138) 2.4 3.0 6.0 0.4 
17 (104,139) 2.9 3.0 5.0 0.2 
18 (104,54) 1.3 1.0 3.9 0.4 
19 (22,104) 1.1 2.0 4.5 0.4 
20 (46,133) 0.7 2.0 4.5 0.4 
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算法 参数 为 : 最 大 迭代 次 数 G=200 , 人 工 狼 个 数 N=20， 表 4 CWPA 最 优 结果 
最 多 使 用 车 辆 数 上 =9 ， 探 狼 比 例 因 子 a=4 ， 最 大 近邻 度 判 Table 4 Optimal consequence of CWPA 
定常 数 c=6 ， 头 狼 变 异 次 数 m=3 ，A=05 。MRGBA 和 车 辆 号 配送 路 径 
PHACA 参数 参见 文献 [3] 和 文献 [11]， 每 种 算法 独立 运行 25 1 0-4-15-9-0 
次 。 2 0-7-11-14-18-0 
图 4 为 CFWA 与 FWA 的 欠 代 优化 对 比 。 可 以 看 出 CFWA 3 0-5-3-0 
在 100 代 时 开始 收 剑 ， 而 FEWA 的 收敛 代数 为 160 代 左 右 ， 4 0-12-19-6-13-0 
并 且 CFWA 的 收敛 效果 要 明显 好 于 FWA。 表 2 记录 了 本 文 5 0-10-20-17-2-0 
法 与 其 他 文献 中 的 两 种 改进 方法 求解 多 目标 VRPTW 的 统 0 
计 结 果 。 在 最 优 车 辆 数 方面 , CWPA 求解 结果 为 6, 优 于 WPA 


表 5 MRGBA 最 优 结 果 
Table 5 ”Optimal consequence of MRGBA 
车 辆 号 配送 路 径 


以 及 文献 [3] 提 出 的 MRGBA 算法 和 文献 [11] 提 出 的 PHACA 
算法 ， 在 最 优 总 距离 方面 ，CWPA 求解 结果 为 1 005.7 km， 
分 别 比 MRGBA 和 PHACA 减少 7.75% 和 13.11%。CWPA 求 


解 的 平均 车 辆 数 以 及 平均 总 距离 也 优 于 WPA、MRGBA 和 2 0-4-15-9-0 
PHACA 求解 结果 ,说明 CWPA 具有 更 好 的 稳定 性 ， 在 平均 2 0-11-14-3-0 
收敛 时 间 方 面 , WPA<CWPA<PHACA<MRGBA, 但 WPA 求解 2 0-18-7-5-0 
精度 较 低 ， 总 体 来 说 CWPA 的 求解 效率 更 高 。 4 0-12-19-6-13-0 
表 2 CWPA 与 WPA、MRGBA、PHACA 求解 结果 对 比 5 0-1-16-8-0 
6 0-10-20-17-2-0 


Table 2 Results comparison among CWPA、 WPA、 MRGBA and 


表 6 PHACA 最 优 结果 


PHACA 
WPA CWPA MRGBAM) PHACAID Table 6 Optimal consequence of PHACA 
最 少 车 辆 数 /veh 8 6 6 7 车 辆 号 配送 路 径 
平均 车 辆 数 /veh 8.36 6.76 6.70 7.36 0-7-3-0 
最 短 总 距离 km 1164.2 1005.7 1090.2 1157.4 2 0-15-4-9-0 
平均 总 距离 km 1290.6 1152.4 1164.1 1187.2 ee 
平均 收敛 时 间 /s 。 51.90 55.17 57.66 56.25 Se 
表 3~6 分 别 为 WPA、CWPA、MRGBA 和 PHACA 各 运 0-10-13-20-6-0 
行 25 次 所 得 的 最 优 结 果 。 图 5 为 本 文 所 提 CWPA 求 得 的 最 和 


优 路 径 图 。 
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Fig.4 Evolution diagrams of CWPA and WPA 怖 攻 
S 图 5 CWPA 最 优 路 径 图 
表 3 WPA 最 优 结果 
Table 3 Optimal consequence of WPA ge pa Ome ol oN 
车 辆 号 配送 路 径 3.3 仿真 实验 2 
pe 为 了 进一步 验证 本 文 所 提 CWPA 算法 求解 不 同类 型 和 
和 a 规模 问题 的 性 能 ， 选 取 Solomon 标准 测试 数据 集中 的 六 个 基 
本 二 条 准 问题 (R1-02、R2-06、C1-04、C2-08、RC1-05、RC2-07) 
Dr 进行 仿真 实验 .CFWA 参数 设置 为 : 最 大 迭代 次 数 G=1000， 
人 人 工 狼 个 数 w=60 ， 最 多 使 用 车 辆 数 依次 为 25、8 、20、8、 
5 0-5-11-0 20、8， 探 狼 比例 因子 a=4 ， 最 大 近邻 度 判 定常 数 cv =6 ， 头 
6 0-10-13-20-6-0 狼 变 异 次 数 m=3 ， 上 =05 。 对 每 个 问题 独立 求解 50 次 ， 
2 0-19-12-0 与 文献 [10] 提 出 的 动态 混合 蚁 群 优化 算法 DHACO 和 文献 [20] 
8 0-8-16-0 提出 的 并 行 禁 忌 搜 索 算法 ITSH 求解 结果 进行 对 比 。 


表 7 CWPA 与 DHACO、ITSH 求解 结果 对 比 
Table7 Results comparison among CWPA、 DHACO and ITSH 
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DHACO0O9 ITSHPOI CWPA 
基准 问题 

OTD OTV MCT/s OTD OTV MCT/s OTD OTV MCT/s 

R1-02 1486 17 68 1487 17 - 1486 17 67 

R2-06 935 3 11 995 3 - 932 3 13 

C1-04 825 10 41 825 10 - 825 10 39 

C2-08 588 3 19 588 3 - 588 3 16 

RC1-05 1579 14 82 1629 13 - 1594 13 77 

RC2-07 1061 3 23 1061 3 - 1061 3 18 
表 7 为 本 文 CWPA 以 及 文献 [10] 提 出 的 动态 混合 蚁 群 优 


化 算法 DHACO 和 文献 [20] 提 出 的 并 行 禁忌 搜索 算法 ITSH 对 
各 算 例 求解 的 最 优 总 距离 (optimal total distance, OTD ) 和 最 
优 车 辆 数 (optimal total vehicle, OTV ) 和 平均 收敛 时 间 (mean 
convergence time, MCT)。 可 以 看 出 ， CWPA 在 求解 基准 问 
题 R1-02 和 RC1-05 时 都 取得 了 较 好 的 解 , 并 且 在 求解 R2-06 
时 取得 了 比 DHACO 和 ITSH 更 优 的 解 ， 在 求解 RC1-05 时 ， 
求 得 的 最 优 路 径 总 距离 大 于 DHACO 所 求 结果 ， 但 均 优 于 
ITSH 所 求 结果 , 在 求解 C1-04 和 C2-08 时 取得 了 与 DHACO 
和 了 ITSH 相同 的 最 优 解 。 
表 7 还 可 以 看 出 CWPA 求 解 大 部 分 算 例 时 的 平均 收敛 
时 间 均 优 于 DHACO ,在 求解 RC2-07 算 例 时 相 比 DHACO 具 
有 明显 的 时 耗 优 势 ， 在 求解 R2-06 算 例 时 的 平均 收敛 时 间 多 
于 DHACO， 但 是 CWPA 求 得 的 最 优 路 径 距 离 更 优 ， 说 明 
CWPA 在 求解 不 同类 型 和 规模 的 VRPTW 问题 时 均 具有 较 好 
的 性 能 。 


4 ”结束 语 


本 文 针 对 多 目标 VRPTW 问题 的 求解 ， 首 先 建立 了 以 最 
小 化 总 距离 和 总 车 辆 数 为 目标 的 数学 模型 。 将 狼 群 算法 嵌入 
文化 算法 的 框架 ， 提 出 了 文化 狼 群 算法 ， 利 用 近邻 度 概 率 矩 
阵 建 立 初 始 狼 群 ， 提 高 了 狼 群 的 优良 性 ， 提 取 当 前 最 优 解 的 
编码 信息 作为 信仰 空间 的 知识 来 指导 底层 狼 群 向 着 有 利 的 方 
向 进化 。 通 过 实验 仿真 验证 了 文化 狼 群 算法 求解 多 目标 
VRPTW 问题 的 有 效 性 ， 进 一 步 与 相关 文献 中 的 其 他 算法 进 
行 比较 ， 表 明文 化 狼 群 算法 具有 更 好 的 稳定 性 。 下 一 步 的 研 
究 目标 为 更 合理 的 设置 文化 狼 群 算法 的 参数 ， 并 且 尝 试 将 其 
用 在 拥有 更 多 目标 个 数 的 多 目标 VRPTW 问题 的 求解 中 。 
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